Frequency shift of optical phonons in doped graphene layers 
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We use first-principles density-functional calculations to determine the frequency shift of the 
Aj-K phonon (Raman D band) in monolayer graphene, as a function of the charge doping. A 
detailed DFT study on the electron-phonon coupling and the phonon line width of E2g-r and A'^-K 
phonons are also performed for graphene multi-layers. Furthermore, we explain the experimentally 
observed '1/ (Number of Layers)' behaviour of the Raman G band position after including the 
dynamic response treated within time dependent perturbation theory. 
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Graphene, a two-dimensional honeycomb lattice of sp^- 
bonded carbon atoms 1], is known to exhibit many re- 
markable properties that are of fundamental interest and 
technological relevance. Being the fundamental building 
block for carbon allotropes of other dimensionality, it can 
be stacked into 3d graphite or rolled into Id nanotube. 
Graphene has created a great interest in the scientific 
community due to its excellent mechanical and electronic 
characteristics, scalability to nanometer sizes [2| and easy 
accessibility to optical probes. Moreover, this atomically- 
thin sheet is thermodynamically stable, continuous on a 
macroscopic scale and exhibits high crystal quality. A 
very intriguing feature of graphene is related to the fact 
that its electron transport is essentially governed by the 
Dirac's (relativistic) equation and can be controlled ex- 
ternally [1, 0, [1, @] ■ Using electric field in an FET geome- 
try it is possible to dope graphene layers by changing the 
carrier concentration in the samples. In particular, its 
phonon spectrum can be modified significantly by tuning 
the applied gate voltage 0, S] ■ 

Phonon dispersions of graphene exhibit two Kohn 
anomalies in the highest optical branches at F and K 
(Ei2g mode - Raman G band and A'l mode - Raman D 
band, respectively). The electron phonon couplings for 
these modes are particularly large and negligible for all 
the other modes at T and K ,9]. In this paper, we cal- 
culate the frequency shift of the long wavelength opti- 
cal phonon (E2g-r mode) in graphene as a function of 
the number of layers. Phonon line width and the elec- 
tron phonon coupling (EPC) for the above two optical 
phonons are also extensively studied in graphene multi- 
layer. Furthermore, for the first time, we compute the 
variation of phonon frequency of the A[ mode at K in a 
graphene monolayer, as a function of doping concentra- 
tion. Calculations are done first using a fully ab-initio ap- 
proach within the standard Born-Oppenheimer approxi- 
mation and then time-dependent perturbation theory is 
used to explore the effect of dynamic response. 

Our ab-initio calculations are performed using den- 
sity functional theory (DFT) [13] in the generalized gra- 
dient approximation. We have used the PWSCF [11 1 



implementation of DFT, with Perdew-Burke-Ernzerhof 
(PBE) 12] for the exchan ge c orrelation functional and 
ultrasoft pseudo potential 13| to represent the interac- 
tion between ionic cores and valence electrons. Kohn- 
Sham wave functions are represented with a plane wave 
basis truncated at an energy cut off of 40 Ry. The two- 
dimensional graphene crystal is simulated using a super- 
cell geometry with a vacuum of 7.5 A in the z-direction 
to ensure negligible interaction between its periodic im- 
ages. The Brillouin zone integration is done on a uniform 
36 X 36 X 1 Monkhorst-Pack [14] grid. An electron smear- 
ing of 0.02 Ry with Fermi-Dirac distribution is used to 
accelerate the slow convergence of self-consistent calcu- 
lation in graphene layers by smoothing out the discon- 
tinuities present in the Fermi distribution at zero tem- 
perature. Structural relaxation is carried out in each 
case to minimize the forces acting on each of the atoms 
using Broyden-Flecher-Goldfarb-Shanno (BFGS) based 
method [15|. Phonon frequencies and electron phonon 
couplings are calculated using Density Functional Per- 
turbation Theory (DFPT) [16| at the level of linear re- 
sponse, which allows the exact (within DFT) computa- 
tions of phonon frequencies at any Brillouin zone point. 
The Fermi-energy shift with doping is simulated by con- 
sidering an excess electronic charge which is compensated 
by a uniformly charged back-ground. 

In graphene, the Fermi surface is reduced to two equiv- 
alent K and K' points in opposite corners of the 2D 
hexagonal Brillouin zone where the valence and conduct- 
ing bands touch each other. This leads to the Dirac cone 
spectrum e^. /^(K-|-k) = ihvpk for the tt* and tt bands, 
where vp is the Fermi velocity of the massless (2-1-1) di- 
mensional Dirac fermions and k is a small vector. Within 
this approximation, at zero temperature (T = K), the 
layer charge concentration and the Fermi energy are re- 
lated via 



cr = sign(eF) 



sign(eF) ij? 



X 10.4 X 10^^ cm-2 (eV)"^ (1) 
where hvp = 5.53 eV-A from our DFT calculation, in 
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excellent agreement with the previous DFT result 
and ep = at the tt bands crossing (known as Dirac 
point). For very small doping concentration (order of 
lO^^cm"^), the linearized bands are a good approxima- 
tion. We have used this model of strongly interacting 
two-dimensional Dirac fermions in few layer- graphene 
as an approximation because the Dirac singularity is the 
topological property of the electronic spectra that should 
be stable towards the weak 3d inter-layer coupling even 
in few layer- graphene. That's why the Fermi velocity 
for the few layers of graphene [18| is taken to be the 
same as for a single layer within a good approximation 
for low doping concentration. To check validity of the 
above approximation in graphene multi-layers, we have 
performed the DFT calculations of the electronic band 
structure for monolayer, bilayer, trilayer and tetralayer 
graphene and our results ( Figure 1 ) confirm that the 
above assumption is justified. 
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FIG. 1: The electronic band structure of graphene layers 
determined with ab-initio DFT calculations. 



mined by the double resonance process [iSSIlll. For 
the typical experimentally used laser excitation of 2.41 
eV, the value of AK is 0.0855 (in units of 27r/a). Our 
results of the frequency shifts for the K point phonon 
{AujII''*''=) and the (K - AK) point phonon (Acj^^^'X^k) 
with charge doping a are shown in Figure 2. In both 
panels, the open circles show the results obtained when 
the lattice spacing is kept fixed at its value correspond- 
ing to the undoped case and the filled stars show the 
results obtained when the lattice spacing is optimized 
for each doping concentration. The optimization leads 
to a lattice expansion with electron doping and a lat- 
tice contraction with hole doping, and gives substantial 
change in the phonon frequency shift. Interestingly, Fig- 
ure 2 shows that the frequency shift of the A'^-K phonon 
on the charge doping in a graphene monolayer is very 
different compared to the F point phonon [17]. 
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FIG. 2; Adiabatic frequency shift of the K- Ai and K - AK 
phonon in graphene monolayer as a function of charge doping 
a, with respect to the zero-doping frequency. 



We first determine the total energy and electronic den- 
sity by solving the self-consistent Kohn-Sham equations 
for the equilibrium crystal geometry where the relative 
atomic positions in the unit cell yield zero forces and 
lattice parameters lead to a zero stress-tensor. Once 
the unperturbed ground state is determined, phonon fre- 
quencies are obtained by using DFPT which calculates 
the linear response of the electrons to a static pertur- 
bation induced by ionic displacements. This approach 
is based on the adiabatic (also referred to as the Born- 
Oppenheimer or static) approximation. At zero-doping, 
our computed E2g-r and A'^-K phonon frequencies in 
a monolayer graphene are Wp"*"^(0) = 1553.3 cm"^ 
and u;|^"**^(0) = 1301.7 cm"!, respectively. Infact, the 
phonon exactly at the K point is not the relevant one for 
the Raman D or 2D bands in graphene; rather what is 
relevant is the phonon at (K — AK) where A K is deter- 



Within the adiabatic approximation, above results are 
generically expected on physical grounds. For the long- 
wavelength r phonon (small wave- vector q), a high den- 
sity of electron doping leads to an effective screening 
of the ion-ion interactions, reducing the elastic coupling 
in the lattice and hence leading to a softening of the 
phonons. On the other hand, in the case of the K-point 
phonon with a larger wave-vector which connects two 
points in the Fermi surface, electrons arc no longer able to 
effectively shield the ion-ion interactions and the vibra- 
tional properties shows a giant and sharp Kohn-anomaly 
- a precipitous phonon softening at K point. With dop- 
ing, the change in the Fermi surface moves that Kohn- 
anomaly away from q = K and stiffens K-point phonon 
(Figure 2). 

So far we have used the standard adiabatic approxima- 
tion which is usually justified in ordinary metals where 
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the electronic energy gap between the ground and the ex- 
cited states is larger than the phonon energy so that the 
phonons respond to a time averaged electron distribution. 
In graphene, however, this approximation is inadequate 
because of its unique massless Dirac-like electron band 
dispersion and comparable electronic and phonon energy 
scale. The complete break-down of the adiabatic approx- 
imation in doped graphene has been recently demon- 
strated, both theoretically and experimentally 0, 
Hence, for single-, double- and few layer- doped graphene, 
we next consider the phonon as a dynamic perturbation 
treated within time-dependent perturbation theory. 

Using such a dynamic approach in the context of DFT, 
the self-energy of a phonon mode v at wave- vector q can 
be written as [22], 



transition and {g\)v = Y^i j l5Ki,KjP/4, the sum is for 
the two degenerate tt bands at ep (see Eqn. 6 and note 
24 of Ref. 11). 

Based on the above approximations with the linear 
electron dispersion near Dirac point, we find that the 
change of Raman G-band phonon energy at zero temper- 
ature can be described as, 
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/(k+q)i ^ /kj" 



£(k+q)i ^kj 



(2) 

where Nk is the number of k-points, the sum is over 
k G Brillouin zone, /kj is the Fermi distribution func- 
tion, 77 is a small real number and <?('i(;+q)i kj 
electron-phonon matrix element. Within DFPT, the elec- 
tron phonon matrix element can be obtained from the 
first order derivative of the self-consistent Kohn-Sham 
potential, with respect to the atomic displacements as: 
5(k+q)*,kj = ((k + q)i\SVKs/Suq^\kj)/y/2Mujq^ where 
Uqu is the amplitude of the displacement of the phonon 
!/ of wave vector q, ujqi, its phonon frequency, M is the 
atomic mass, Vks the Kohn-Sham potential and |kj) is 
a Bloch eigenstatc with wave-vector k, band index j and 
energy e^j. 

The dynamic correction to the energy shift arising from 
the electron-phonon interaction is given by the real part 
of the phonon self-energy as: 
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/(k+q)i - /kj 



.e(k-i-q)i - Ckj 



■ CO. 



(3) 

where V takes the principal part of its argument. 

The phonon line width (FWHM) is twice the imaginary 
part of S'i/(q, Wqi,) and can also be obtained from the 
Fermi golden rule [2^: 



^l'' = l^Z^lf(k+q 



{fkj - /(k+q)i) '5(e(k+q)i - Ckj - t^qi.) (4) 

Ideally, in the dynamic case, lu should be determined 
self-consistently. However, considering the dynamic and 
doping effects as perturbation, at the lowest order, we 
can use the adiabatic undoped phonon frequency ujq"'*^'^ 
in place of u. For q near F-point, we can also con- 
sider the electron phonon coupling as |(7(K+k)j,(K+k)iP = 
(.gp)F[l ± 008(26*)], where 9 is the angle between the 
phonon-polarization q and k, the sign ± depends on the 
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where ar 



(5) 

Aq = 5.23 is the equilibrium 
unit cell area. In this case, the energy shift diverges 
logarithmically when the magnitude of the Fermi energy 
becomes the half of the phonon energy and increases in 
proportion to the Fermi energy for \ep\ > ?i(jjQ*°*"^/2 as 
long as the effect of the charge doping can be considered 
as a perturbation. Beyond this, one can easily notice a 
considerable change in ar and cjp "'^'^ with charge doping. 

In order to calculate the dynamic frequency shift of 
Raman G band for graphene layers, we need ar for 
different layers. This is calculated using ab-initio den- 
sity functional theory and results are given in Table I. 
Further, we find that the static zero-doping frequencies 



= (0) 



1301.7 cm~ 



and wf*°"=(0) 



1553.3 cm 



do not change with the number of graphene layers, n. We 
plot the dynamic frequency shift of the Raman G band 
as a function of the charge doping for different graphene 
layers in Figure 3(a). One can readily notice (from Fig- 
ure 2) that for small doping (around 0.4 x 10"^"^ cm~^ ) 
expected in zero-biased unintentionally doped graphene 
on a substrate, the static frequency shift is very small 
and can be treated as a correction ( < 1 cm~^ ) with 
respect to the dynamic one. For such a low doping con- 
centration, the dynamic shift shows a linear dependence 
on 1/n (Figure 3(b)). Interestingly, this is very similar 
to the the experimentally observed [13] dependence of 
the the Raman G band position on the number of layers. 
We, therefore, suggest that the frequency shift of the Ra- 
man G mode on the number of layers is dominantly due 
to the unnitentional charge doping which can be under- 
stood based on the dynamic effects. 

Moreover, our DFT calculations (Table I) show that 
the phonon linewidth changes very little with the num- 
ber of graphene layers. This result is expected on phys- 
ical grounds. In graphene layers, the phonon linewidth 
due to the electron phonon coupling obtained from the 
imaginary part of the phonon self-energy (or equivalently 
from the Fermi Golden rule) is 0, [l3| 
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where n is the number of graphene layers. We find that 
(Q!r)„ n is almost constant with the number of layers. 
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Number of a l[T = GK] {cm~'^) 

Layers, n Qr Qk 7r 7k 

1 4.441 X 10"'' 13.33 x 10^'' 11.70 SlO 

2 2.223 X 10"^ 6.68 x 10"^ 11.56 8.04 

3 1.457 X 10*^ 11.43 

4 1.085 X 10"^ 11.32 

5 0.889 X 10~^ 11.23 



TABLE I: Electron-phonon coupling strength, a and phonon- 
line widths, 7 (cm~^) at T and K for graphene layers. The 
values are computed using density functional theory, a de- 
creases very fast with increasing number of layers while the 
change in 7 is very small. 

leading to a very small variation of phonon line width 
with number of layers. At zero temperature, we get 
7(T=0K) = ^("r)„n?»^g'°"- e(fta;^*"*^^-2|eF|). It implies 
that the broadening is nonzero only for IefI < fiu;g*°**'^/2 
when the scattering process is allowed by the Pauli ex- 
clusion principle. 



will be reported elsewhere very soon. 

In conclusion, we have computed, from first-principles, 
the frequency shift of the A'j^-K phonon (Raman D band) 
in monolayer graphene, as a function of the charge dop- 
ing. A detailed DFT study on the electron-phonon cou- 
pling and the phonon line width of E2g-r and A'j^-K 
phonons have also been performed for graphene multi- 
layers. Furthermore, we explain the experimentally ob- 
served 'l/(Numbcr of Layers)' behaviour of the Raman 
G band position after including the dynamic response 
treated within the time dependent perturbation theory. 

S.K.S. gratefully acknowledges the financial sup- 
port of Research Fellowship from Council of Scientific 
and Industrial Research (India) under Award Number 
9/79(913)/2002-EMR-I. A.K.S. thanks Department of 
Science and Technology (India) for financial support. We 
also thank Prof. Prabal K Maiti for providing access to 
his computer resources. 




Electron concentration ( 10 cm ^ ) 




^ 1 / Number of Graphene Layers, n 

FIG. 3: (Color online) (a) Dynamic frequency shift with 
charge doping in graphene multi-layers and (b) Frequency 
shift of Raman G band as a function of the number of layers. 

Finally, we discuss the effect of the charge doping on 
the K-point phonon. Recent experimental observation [8] 
indicates that although the width and the peak position 
of the Raman G and D* bands ( second order of K point 
phonons ) exhibit similar doping dependence, the magni- 
tudes of the changes are only around 10 percent of the G 
band. This small dependence suggests that for the large 
wave vector phonon (K point phonon) there is a cancel- 
lation of the large inter-band and intra-band contribu- 
tion [25] in the dynamic response and the static DFT 
frequency shift (see Figure 2) plays the major role. We 
are also studying the layer dependence of the K-point 
frequency shift in doped graphene multi-layers and that 
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